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Abstract 

In low energy SUSY theories exchange of gluino and squark with left-right 
mixing can produce a large gluonic dipole interaction. In this paper we study 
the effects of this interaction on B ^ Kir using QCD improved factoriza- 
tion method. The Standard Model predicts a smaller branching ratio for 
— > K^TT^ than experimental measured one. We find that within the pa- 
rameter space allowed from B jXs constraint, the SUSY dipole interac- 
tion can enhance this branching ratio to agree with the experimental mea- 
surement. Combining recent data for all the four B^ K^ir^ , K^-k^ and 
B~ — > K~TT^ , K^TT~ decay modes, we find that the allowed parameter space 
is reduced significantly compared with that using B Xg'j data alone. It 
is found that even with these constraints, the predictions for CP violation in 
these modes can be dramatically different from those of the SM predictions. 
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I. INTRODUCTION 



There have been considerable experimental and theoretical efforts to understand the 
properties of B systems. These studies have provided important information about B decays 
and CP violation. At the quark level the relevant Hamiltonian for B decays in the Standard 
Model (SM) is well understood. When going beyond the SM, there are new contributions. 
These new contributions can modify or even completely change the SM predictions In 
SUSY theories with flavor changing interaction in the squark sector, it is possible to generate 
large effects on hadronic B decays while their effects on other processes are small 0,^. In 
particular exchanges of gluino and squark with left-right mixing can enhance the gluonic 
dipole interactions of the forms qa^^ G^''^ {l±^^)h by a large ratio factor of gluino mass nig to b 
quark mass m^, mg/mi,, compared with the SM prediction. Due to this enhancement factor, 
even a tiny coefficient for the associated flavor changing squark-gluino-quark interaction, a 
large gluonic dipole interaction can be generated. 

A large gluonic dipole interaction can affect B decays significantly. It may help to explain 
the large branching ratios observed for B XsT]', although theoretical understanding is poor 
[^J^. It can also change theoretical predictions for other charmless hadronic B decays, such 
a.s B ^ K7i,7nT,(f)K P,^. Using the recently measured branching ratios for B Kn^Tnr, 
one may be able to constraint the allowed parameter space which can generate a large gluonic 
dipole interactions and to provide interesting information about models beyond the SM. In 
SUSY models the same left-right squark mixing parameters can also generate a photonic 
dipole interaction which can induce B Xg'y and B Xd'~f. At present B — > Xg'y has 
been observed, but not B X^'y. 

In this paper we study in detail constraints on SUSY gluonic dipole interaction using 
data from B^ — > K^n^, K^tt^ and B^ K^n^, K^n^ and also the measured branching 
ratio oi B Xg'j. Previous studies of gluonic dipole interactions on B decays were based 
on naive factorization approximation. Here we will use QCD improved factorization method 
I^J^ which improves the analysis on several aspects, such as the number of color, the gluon 
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virtuality, the renormalization scale, scheme dependencies and etc. We find that the recent 
data on B Ktt decays can reduce, significantly, the allowed parameter space compared 
with using B —>■ Xs'~f data alone, while still allow large deviations from the SM predictions 
for branching ratios and CP asymmetries in these decay modes. 

The paper is arranged as follows. In section II, we discuss gluonic dipole interactions 
in low energy SUSY model with left-right squark mixing, and obtain the decay amplitudes 
for B Kit decays using QCD improved factorization. In section III, we first update con- 
straints on the SUSY gluonic dipole interactions from B — > Xg'y, and then study constraints 
from B — > Kit decays. In section IV, we study CP rate asymmetry in B —>■ Kir. And in 
section V, we draw our conclusions. 

II. SUSY GLUONIC DIPOLE CONTRIBUTIONS TO B ^ Kn 

In the SM the effective Hamiltonian for charmless photonic and hadronic B decays with 
AS* = 1 at the quark level is given by 

H,ff = ^ lvubV:XciO, + C2O2 + E c,0,) + V,tK: (^^0^ . (1) 

Here Oj are quark, gluon and photon operators and are given by 

Oi = {SiUi)v-A{Ujbj)v-A, O2 = {SiUj)v~A{Uibj)v-A, 

03(5) = {sibi)v-A'Y{qjqj)v-{+)A, 04(6) = {sibj)v-A'Y{qjq'i)v-i+)A, 

q' q' 
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q' 



On = ^m,-s,a^''GlXK^ + O12 = -^m.s^a^'' F,,{1 + 75)6., (2) 

where {V A){y ± A) = ± 75)7^(1 ± 75), q' = u, d, s, c, b, Cqi is the electric charge 
number of the q' quark, is the color SU(3) generator normalized as Tr{T°'T^) = 6"'^/2, i 
and j are color indices, and Gfj,u and F^^ are the gluon and photon field strengths. 

The Wilson coefficients q have been calculated in different schemes ||^. In this paper we 
will use consistently the NDR scheme. The values of q at ^ nib with QCD corrections at 



NLO are given by [0 

ci = 1.081, C2 = -0.190, C3 = 0.014, a = -0.036, cg = 0.009, ce = -0.042, 

Cj = -O.Ollaem, Cs = O.OeOOfem, Cg = -1.254aem, Ciq = 0.223aem. (3) 

And at LO are given by [0 

ci = 1.117, C2 = -0.268, C3 = 0.012, C4 = -0.027, C5 = 0.008, Cg = -0.034, 

C7 = -O.OOlaem, Cs = 0.029aem, Cg = -1.276aem, Clo = 0.288aer„, 

cii = -0.151, C12 = -0.318. (4) 

Here aem = 1/129 is the electromagnetic fine structure constant. We use LO Wilson coeffi- 
cients with the matrix elements enter the decay amplitudes at next-to-leading order. 

In SUSY models, exchanges of gluino and squark with left-right squark mixing, can 
generate a large contribution to Cii^i2 at one loop level. In general these contributions can 
generate a gluonic dipole interaction with the same chirality as the SM one, as well as with 
an opposite chirality as the SM one, that is, an interaction similar to On, 12 but with 1 -|- 75 
replaced by 1 — 75. It is difficult to carry out an analysis in the full parameter space. We 
will first consider the new contributions with the same chirality as the SM one setting the 
opposite chirality one to zero, and then the opposite case. 

The effective Wilson coefficient for clYu obtained in the mass insertion approximation 
is given by, for the case with the same chirality as the SM ones [0], 

C12 \"^g) - 77—2 Y7-T7^—^^y^mh 



X 

3(1 



Ga{x) = —-^-—[22 - 20x- 2x2 + 16xln(x) - x^ ln(x) + 9 ln(x)], 

- X) 

Ax 

Fq{x) = -— y[l + Ax - 5x^ + 4:xln{x) + 2x'^ln{x)], (5) 

where 5ff^ is the mixing parameter of left and right squarks, Xgq = m^/m?^ is the ratio of 
gluino and squark mass. 



The coefficients cl\^i2{mg) at are given by [|T0 



susy / \ susy / \ susy / \ 2 susy / \ \ t 2\ susy / \ /c\ 

where r] = {as{mg) / as{mt)y^'^^ {as{mt) / as{mb)Y^'^^ . 

In the SM, cu^u are proportional to mb/rn^r. From the above expressions, it is clear that 
the SUSY contributions are proportional to l/rrig. If nig is of order a few hundred GeV, 
there is an enhancement factor of rrig / mh{ni^ / m"^) for the SUSY gluonic dipole interaction. 
In this case even a small 5ffj:j can have a large effect on B decays. 

To obtain clYu for opposite chirality case, one just adds in two more operators similar to 
Oil, 12 but with 1 + 75 replaced by 1 — 75 and the Left-Right mixing parameter S'^j^ replaced 
by the Right-Left mixing parameter (5jjj^. 

We follow Ref. to obtain the B Kit decay amplitudes. They are given by 

+VctK:{-al + ^al, - RK{al - \al)] - 9z,f^fj^iy^^v:, + K,K:)(63 - 



+K.V;:(a^ + RK{al + al + a^J + a^,)] + ^^f^fMV^,v:,{h + h + bf"^) + K.Ktl^s + bf"^)], 
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Here Rk — 2m\/msmb. ai and bi coefficients are related to tlie Wilson coefficients. In 
the above we have neglected small contributions from Oi2- Including the lowest as order 
corrections, af are given by 

C2 Cpiras 



02 



04 = C4 + 



C3 



p 



ay 



. C5 

^« + iV 



1 + 

1 + %^v. 

47r 



+ 



-H, 



Ktv, 



+ — — -HttK, 



1-6 



+ 



+ 



6 ■ 



CpOts 



CfCXs Pk,3 

An ' 



+ 



H 



Ktt, 



An 



+ 



Q^em ^K,3 

9n Nr ' 



1 + 



Oio = Cio + — 



1 + 



47r 
An 



Cio CFTTtts 

+ ^r;^ — — rl-n-K, 



K 



+ 



r>P,EW 

9n Nr 



+ — TT — ri 



Nc N, 



where Cp = {N^ - l)/(2A^c), and N^ = 3. The quantities H^.m,, ^k,2, ^x,3, PrS 
and -P^'f ^ are hadronic parameters that contain all nonperturbative dynamics. 
The vertex corrections V^''(M=7r,K) are give by 



Ktt: 



(8) 



yP,EW 



Vm 
9{x) 



121n— -18+/ dxg(xUMix), 
jj, Jo 

12 In— -18+ C dxgil- x)(Pm{x), 
II Jo 



= 3- 



/i 

2x 



1 — X 



Inx — Sin 



+ 



2Li2(x) -ln^x + 



2 In a; 
1-x 



- (3 + 2in) -{x ^1-x) 



(9) 



where Li2(a;) is the dilogarithm. (j)M{x) is the leading- twist light cone distribution amplitude. 
This distribution amplitude can be expansion in Gegenbauer ploynomials. We truncate this 
expansion at n = 2. 



3m I 



x,ii) = Qx{l -x) 1 + af (/x)Cr''(2a; - 1) + <(/^)C^' 



(3/2), 



Ml 



.(3/2)(2x-l) 



(10) 
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where Ci'^'^\u 



3m and Cf'^\ 



u) = |(5u —1)- The distribution amphtude parameters 
afl^ for M = K,7r are: af = 0.3, af = 0.1, < = 1 and = 0.1. 

The penguin contributions -P|- 2, Pr^^ given by 

2 



^^,2 = ci 



with 



3ln- + --G;.(.,J 



+ C3 



^ln^ + ^-GK(0)-G,,(l) 



+ (c4 + ce) 



1^ In - (n^ - 2)G;,(0) - Gk{s:) - Gk{1) 

3 yU 



dx , , . 

-2cii / (Pk[x) 

Jo I — X 

(Ci + N^C2) 



2ci2 / -^^(pxix) 
1 — a; 



+ (c4 + ce) 



4: rrib 2 a . . 
-In — + - - Gk{Sp} 
3 /i 3 



+ C3 



^ln^ + |-G/.(0)-G,.(l) 
3 /i 3 



1^ In !^ - (n^ - 2)G^(0) - - Gk{1) 



2cii, 



- m GK[Sr,) 

3 2 ^ 



3c 



12, 



/ dxGis — ze, 1 — 

JO 



dxG{s — ie,l — x) 



G{s, x) = —4 / duu{l — u) ln[s — ^(1 — u)x] 
Jo 



(11) 



(12) 



where 4>p{x) = 1. 

The hard-scattering contributions H^^k, H'^j^, and H^-k contain many poorly know pa- 
rameters,we follow M use H^^k = 0.99 and 



Hir^ = H, 



(13) 



where 



R 



ttK 



^0 



(0)/. 



Fo^-^{^)fK 



(14) 



All scale-dependent quantities for hard spectator contributions are evaluated at = y/Ahrub 
with A;, = O.SGeV. 



The annihilation coefficients are 



c 

&r=^M+c8A^]. (15) 

c 



with 



'^A - 4 + y ] + 2rlXl 



A\ = A\ = Tcas 

Ai = l27ia,rl{2Xl-XA), (16) 



where = 2^yjm},, Xa = Jody/y parameterizes the divergent endpoint integrals, We 
use Xa = ln(mB/A/j). All scale dependent quantities for annihilation contributions are 
evaluated at fih- 

In the numerical evaluation of these expressions we follow Ref. [0] consistently drop 
high-order terms in the products of the Wilson coefficients with the next-to-leading order 
corrections. I.e., we evaluate all the corrections beyond naive factorization with leading 
order Wilson coefficients. 

For the case where the SUSY contributions have the same chirality as the SM one, 
cii = cfi^ + cli^^ . For the opposite chirality case, one needs to change cn to cfl'-' — cli^^ due 
to the replacement of 1 + 75 by 1 — 75. 



III. CONSTRAINTS ON Cn FROM B 7X5 AND B ^ Kn 

A. Constraints from B Xg'-f 

We now study the constraints on gluonic dipole interaction using B Xs'y. To this end 
we define = |ci2*^(mb)/cf2^(mb)| and Vg = |ci"*^(mb)/cf/^'^(mb)|. cff [2(^11,) denote SM 
contributions, which are given in the Eq.(||). Using Eq.(|D and replacing cli^i2 with Eq.(|), 
we obtain 
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[v 



Go(x, 



gq, 



6 Cii 



(17) 



r^/rg is not sensitive to gluino mass, but strongly depends on Xgq. In Fig.(|l|) we show r^/r^ 
as a function of Xgq for different gluino masses. Note that both and have a common 
CP violating phase 5 which is equal to the phase difference of the SM phase and the SUSY 
phase of 5f^j^. 
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FIG. 1. Ratio of rg/r^ vs. Xgg. The solid, dashed and dotted lines correspond to gluino mass 
of 1000 GeV, 500 GeV and 200 GeV, respectively. 

For a given r^, Vg is known as a function of Xgq. Therefore constraints on from 
B Xgj can be translated into constraint on Vg. The constraints on is given by 



Br{B X,7) = Br{B Xeve) 



(18) 



I Kb 1 2 ng{mc/mb)r] 

where c[{^ {M^) = cff{mb){l + r^e'^), g{z) = 1 - 8^^ + 8z^ - - 24:zHnz, and r/ = 
1 — 2/(r, 0, 0)as(mb)/37r with /(r, 0,0) = 2.41 [|ll|]. Here we have used the leading order 
result. The use of next-leading order result ||12[ will slightly change the results (at the level 
of 10%), but will not change the main conclusions. 



For the opposite chirality case, one needs to replace |c 



e//|2 
12 I 



i5\\2 



in the 



above by |ci2'^p 



.SM\2i 
12 I I 



1 + r^) because the SM and SUSY contributions have different 
chiralities. Note that in this case r^^g also have a common phase equal to the phase difference 
of SM and SUSY due to which will be also indicated by 6. 
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In the numerical analysis we will use Br{B Xeve) = 10.4%, rrih = 4:.2GeV, rric 



l.SGeV, and l^^^^feP/l Vc^p = 0.95. Using the experimental allowed range Br{B Xg'j) = 
(2.96±0.35) X 10"^ averaged from CLEO, ALEPH, and Belle g, one can obtain the allowed 
region for by Eqs.(13) and (14). The 95% CL. region of r-y is shown in Fig.(0). Combining 
information from Figs.(0) and we obtain the allowed region for Vg in Fig.(|^). We see 
that cl\^^ can be considerably larger than the SM. We note that the bounds of and Vg in 
the opposite chirality case can be up twice depend on the choose of the parameters mc/mb 
and Ci2. It is, however, not change our results too much in the later discussion. 

2.5r 




50 100 150 200 250 300 350 
5 

FIG. 2. Constraint on from B Xg'j. For the case with the same chirality as the SM one, 
the 95% C.L. allowed region is indicated by the shaded region. For the opposite chirality case, the 
% C.L. allowed region is below the dashed line. 



B. Constraints from B Ktt 

It is clear from discussions in the previous section B Xs'~f can constrain clY^ , but 
still allow large deviations from the SM prediction. Within the allowed regions for Vg, rare 
B — s> Ktt decays can be dramatically different from the SM predictions. Therefore using 



experimental data on B ^ Ktc given in Table I |Tj], one can further constrain the allowed 
regions for r^. 

The branching ratio for B —>■ Ktt with SUSY contribution to cn can be easily obtained 
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FIG. 3. Constraint of rg in 95%C.L., with 200 < m~ < 1000, 0.1 < Xgq < 10. The shaded 
region and the region below the dashed hne are the ahowed regions with the 95% C.L. for the case 
with the same and opposite chirahties, respectively. 

by using cn = cff^{l + r^e"') in Eq.(8) for the SUSY contributions with the same chirahty 
as the SM one, and by using cn = cf/^^(l — r^e"') with opposite chirality case. To finally 
obtain the branching ratios, we need to know the form factors, F^^^, the quark masses, 
the parameters of distribution amplitudes, the CKM parameters Vub, Va, and the phase 7. 
There are several theoretical calculations for the form factors, we will use Fq^'^ = 0.28 
and Rt^k = 0.9 given in Ref. [0. The quark masses mh{mb) = 4:.2GeV, rricirni) = l.SGeV, 
m,{2GeV) = llOMeV and (m„ + md){2GeV) = 9.1MeV will be used for illustration. For 
the magnitudes of the CKM parameters we will use \Vcb\ = 0.0402 and iKit/Kfel = 0.090 and 
treat 7 as a free parameter. 

To see how the SUSY gluonic dipole interactions can change the SM predictions in detail, 
we first study some special cases for On with the same chirality as the SM one, and then 
consider the combined constraints. For the special cases we take determined using the 
central value for B —>■ Xg'-f, and take the corresponding Vg with Xgg = 10 to minimize the 
effects. For the CP violating phase we consider three scenarios: a) r^e**^ is real and the phase 
7 is the only CP violating phase; b) the phase 7 is set to zero and the new contribution has 
a phase 6 which can vary from to 2tt; And c) the phase 7 is fixed at the current best fit 



value 7 = 66° |jT6| in the SM and let the phase 6 to vary from to 27r. 
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Branching ratio 


data 


Average 


Br{B K+TT-) 


CLEO 

Belle 

BaBar 


17.2l|;^ ± 1.2 

iQ Q+3. 4+1. 5 

— o.z — U.O 

16.7 ± 1.6 ± 1.3 


17.3 ± 1.5 


BR{B K-7r°) 


CLEO 

Belle 

BaBar 


11 R+3. 0+1.4 
il.U_2 7_i 3 

16.3+11+11 

~o.o 1 .o 

w.8tl:l ± 1.0 


12.1 ± 1.7 


Br{B K^TT-) 


CLEO 

Belle 

BaBar 


18.21^1; ± 1.6 

1 Q 7+5.7+1.9 
-■■"J- ' -4.8-1.8 

18.2l||j ± 2.0 


17.3 ±2.7 


Br{B K^TT^) 


CLEO 

Belle 

BaBar 


14 fi+5.9+2.4 
J-^-0-5. 1-3.3 

ifi n+^-2+2-5 

8.2^2,7 ± 1.2 


10.4 ±2.7 



TABLE I. Experimental data of -B ^ Kir decays from CLEO, Belle and BaBar [14 

The results are shown in Fig.(^). For comparison, we first show the SM predictions for 
the four B Kn decays in Fig.(§.i). The 7 phase seem to be allowed only in the small 



regions around 100° ~ 109° and 251° ~ 260°, which are different from the globe fitting [|16 
Of course this situation can be different when different parameters are used 0. we will take 
the value obtained here for illustration. 

The situation can be dramatically changed if SUSY gluonic dipole interactions are in- 
cluded. In case a), there are two regions of solutions of VgC^^, one is small negative and the 
other is negative with large magnitude. The small negative solution r^e*'' is too small to 
influence B — > Kit decay. For the large negative r^e**^ we find that the gluino contributions 
is 5.6 times larger then the SM one, but have opposite sign to the SM one. In this case, all 
B Kn decays are too large as can be seen from Fig.(^.ii). If a smaller Xgg is used, the 
situation is worse. SUSY contribution with no phase {6 = 0) is not favored. 
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In case b), when the phase S is small, the branching ratios are almost the same as those 
in the SM which increase with the phase 6 to maximal at 5 = 180°. There is a large gap 



between B 



TX 



and B K+^- shown in Fig.(|.iii). B K^ti^ and B ^ R-n^ 



are overlap. The regions with branching ratios for B Kit to be within experimental 2a 
ranges are located between 117° ~ 120° and 240° ~ 243°. 

In case c), the four branching ratios are shown in Fig.(^.iv). In this case it is, again, 
possible to make all decays into experimental 2cr ranges. The overlap regions for 6 are 
located between 109° ~ 121° and 243° ~ 261°. 



{\)BR{B Ktt) vs. 7, SM prediction. 
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{u)BR{B Kir) vs. 7, is real. 
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(iii) BR{B Ktt) vs. 5, 7 = 0° 



(iv) BR{B Ktt) vs. 6,-/ = 66°. 
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FIG. 4. Branching ratios oi B ^ Ktt. The line segments with shade indicate the branching 
ratios which are within the experimental 2a regions. 

From the above discussions, we see that SUSY gluonic dipole interaction can improve 
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agreement of theoretical predictions and experimental data. Both CP violating phases 7 
and 6 can affect the branching ratios significantly. 

Since the gluonic dipole interaction can have significant influence on B Kir branching 
ratios, B — > Ktt decays can, therefore, also be used to constrain SUSY gluonic dipole 
interaction. To obtain the restricted regions, we scan the allowed region of Vg obtained earlier 
from B Xs'j and the SUSY parameter space 0.1 < Xgq < 10, 200 GeV < nig < 1000 GeV. 
The allowed regions for Br{B Kn) within 2a of experimental values are shown in Fig.(^. 
In obtaining the allowed regions in Fig.(^), we treated 7 as a free parameter varying in the 
range 0° ~ 360°. Figs.(^a-d) show the allowed regions from each B Ktc decay. All four 
decays constrain Vg to be less than 7.2. B^ K^n^ provides the most powerful constraint 
on allowed region for r^. The regions with > 4.4 are ruled out at 95% C.L.. 

In Fig.(|^.e) we use all four B Kit decays to constrain Vg with 7 varying from 0° to 
360°. Tg is constrained to be less than 4.4. 6 between 0° ~ 80° and 280° ~ 360° are allowed 
with small Vg disfavored due to the reason that it make the branching ratio of B ^ Kn 
smaller. We also show the constraint with 42° < 7 < 87°, which is the 95% CL. of 7 from 
the fit of unitarity triangle in the SM |T^ in Fig.(^f). In this case we see that the allowed 
regions of Vg and 6 are further reduced, 0.3 < Vg < 4.3 and 6 is between 100° ~ 274°. 
Non-zero Vg and 6 fit data better than SM. 

The constraints on the SUSY gluonic dipole interactions with opposite chirality to the 
SM one can be easily obtained by using cn = cff''^{l — r^e*"^) in Eq.(8). The combined 
allowed regions on Vg, under the same conditions as for the previous discussions, are shown 
in Fig.(^. Data from B Kn, again, can further constrain the allowed parameter space 
compared with constraint from B Xg'y alone. The region between 100° < 6 < 260° is not 
favored. 
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(e) All B ^ K-K mode, 0° < 7 < 360° (f) All B ^ Ktt mode, 42° < 7 < 87' 




5 5 



FIG. 5. Constraints on rg using B Xgj and B Ktt decays. The solid lines are the 
boundaries of the constraint from B Xgj and the dotted regions are the allowed regions by data 
from B Kit at 2a level. 
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FIG. 6. The allowed regions for SUSY gluonic dipole interaction with opposite chirality as 
the SM one. The dotted regions are the allowed regions from B Ktt constraints, and the region 
below the dashed line is the upper bound from B — > constraint. 

IV. CP VIOLATION IN B ^ Ktt WITH SUSY CONTRIBUTIONS 

From previous discussions, it is clear that SUSY gluonic dipole interactions can affect 
the branching ratios of i? ^ Ktt significantly. It is interesting to note that large new CP 
violating phase 6 is still allowed which may affect CP violation in these decays. In this 
section we study CP violating rate asymmetry 



with SUSY gluonic dipole interactions. 

The allowed CP asymmetry Aasy for each of the B Kit decay is obtained by using 
the allowed regions of parameter space constrained from the previous section. The results 
for SUSY gluonic dipole interactions with the same chirality as the SM one and opposite 
one are shown in Fig. (0) and Fig. (P), respectively. The reference SM predictions as a 
function of 7 are shown as solid lines in Fig.(|^). When SUSY gluonic dipole contributions 
are included, the asymmetries can be dramatically different because Vg and 6 both can be 
large. The scattered dots in Fig. (|^ above and bellow the SM predictions correspond to 
the regions on 6 > 180° and 6 < 180° in Fig.(|^), respectively. We see clearly that the 
predictions can be much larger than the SM predictions. For example, with Vg = 3, 6 = 120° 



F-r 



(19) 



asy — 



F + F' 



16 



and 7 = 66°, the branching ratios for K^n , K vr*^, K^n , and i^^vr", are predicted to be 
(16.8, 10.7, 21.8, 9.2) x 10^^ which are within the 2a region of data, and the asymmetries 
are predicted to be -0.04, 0.00, -0.06, -0.10, respectively. With = 3, 5 = 250° and 
7 = 66°, the branching ratios are (19.0, 12.1, 20.0, 8.1) x 10~^, and the asymmetries are 
0.11, 0.12, 0.09, 0.07, respectively. 

For SUSY gluonic dipole interaction with the opposite chirality, the allowed CP asym- 
metry Aasy for each of the B — » JCvr decay are similar as the above examples. The results 
are shown in Fig.(P). 

A particular interesting case is that for B~ — > K^tt^. The SM predicts that \Asym{B~ — 
K^7r^)\ < 1%. In the [0, considering all kinds of uncertainty, \Asym{B^ K^n^)\ is still 
not larger than 2%. With SUSY dipole interaction, \Aasy{B^ K^n^)] can be as large as 
10%. Observation of CP violation in B^ K^tt~ significantly large than SM prediction is 
an indication of new physics. 

V. CONCLUSIONS 

In this paper we have studied the contributions of SUSY gluonic dipole interaction to 
B Kit decays. We found that SUSY gluonic dipole interactions can affect B — > Kit sig- 
nificantly. QCD improved factorization calculation is lower than the current experimental 
branching ratio for B^ —>■ K^tt^. If experimental data will be further confirmed, this is an in- 
dication of new physics. We found that SUSY dipole interactions can improve the situation. 
All four measured B —>■ Ktt decays can be in agreement with theoretical calculations with 
SUSY gluonic dipole interactions and at the same time satisfy constraint from B Xg'y. 

Present data from B Kit can also further constrain the allowed parameter space 
for SUSY gluonic dipole interactions. A large portion of the parameter space allowed by 
B —>■ Xs7 are excluded hj B ^ Kn data. SUSY gluonic dipole interaction coefficients Vg 
and 6 are constrained to be into two narrow regions for both types of dipole chiralities. 

Constraints from B Xg'y and B Ktt still allow large new CP violating phases in the 
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FIG. 7. The CP asymmetry Aasy for B Ktt decays with SUSY gluonic dipole interaction 
having the same chirahty as the SM one. The soUd curves are SM predictions and the dotted 
regions are predictions with SUSY gluonic dipole interaction. 
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(a) Asymmetry of K~^7r 



(b) Asymmetry of i? K 




50 100 150 200 250 300 350 
Y 




50 100 150 200 250 300 350 
Y 



(c) Asymmetry of -B — > K^tt 
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FIG. 8. The CP asymmetry Aasy for B Ktt decays with SUSY gluonic dipole interaction 
having the opposite chirahty as the SM one. The sohd curves are SM predictions and the dotted 
regions are predictions with SUSY gluonic dipole interactions. 
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SUSY gluonic dipole interactions. This allows very different predictions for CP asymmetries 
for B —>■ Kn decays. In particular, with SUSY gluonic dipole interactions the value of CP 
asymmetry for K^n^ can be as large as 10% which is much larger than the SM 

prediction of less than 1%. This can provide an important test of new physics beyond the 
SM. CP violation in other modes can also be larger than the SM predictions. 

Finally we would like to comment that with SUSY gluonic dipole contributions, some re- 
lations predicted in the SM may be violated. One such example is the rate difference, defined 
as A = f -F. In the Standard Model, one has A{B° -> K-n+) = -{fk/UfA{B^ n+n-) 



17|,if kaon and tt have the same wave function distribution amplitude. With different dis- 



tribution amplitudes the relation is violated at 10% level. In obtaining this an important 
property of the CKM matrix element Irn(yubV*^V^lVts) = -Im{VudV*^V^^Vtd) has been used. 
Since the SUSY contributions are proportional to 5^^^^ for B Kit decays, and propor- 
tional to S'j^pi^ for B ^ -KIT decays. In general 6'}^^^ and are not related, the con- 
tributions from SUSY will break the relation A(5° ^ K-n+) = -{fK/f^?A{B^ tt+tt-). 
We have checked with numerical calculations and found that indeed when SUSY contri- 
butions are included there are regions of parameter space where the relation mentioned is 
badly violated. Experimental measurements of these rate differences can also serve as tests 
of new physics beyond the Standard Model. 
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